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Abstract 

To investigate the effects of changes in the latitudinal temperature gradient and the global mean 
temperature on dust concentration in the Northern Hemisphere, experiments with the GISS 
GCM are performed. The dust concentration over Greenland is calculated from sources in 
central and eastern Asia, which are integrated on-line in the model. The results show that an 
increase in the latitudinal temperature gradient increases both the Asian dust source strength 
and the concentration over Greenland. The source increase is the result of increased surface 
winds, and to a minor extent, the increase in Greenland dust is also associated with increased 
northward transport. Cooling the climate in addition to this increased gradient leads to a 
decrease in precipitation scavenging, which helps produce a further (slight) increase in Greenland 
dust in this experiment. Reducing the latitudinal gradient reduces the surface wind and hence 
the dust source, with a subsequent reduction in Greenland dust concentrations. Warming the 
climate in addition to this reduced gradient leads to a further reduction in Greenland dust due 
to enhanced precipitation scavenging. These results can be used to evaluate the relationship of 
Greenland ice core temperature changes to changes in the latitudinal and global temperatures 
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dient and the changes in the average temperature 
may cause a change in the dust signal in Greenland 
by changes in wind speeds, precipitation, or changes 
in circulation patterns. In an attempt to distin- 
guish these alternatives, in this paper we investigate 
the consequences of such changes on dust concentra- 
tions at Greenland. This follows an investigation by 
Rind [1998], who evaluated the effects of a change 
in the latitudinal sea surface temperature (SST) gra- 
dient using the Goddard Institute for Space Studies 
(GISS) GCM. (Here we do not attempt to use a GCM 
to simulate all of the possible processes involved in 
creating the strong dust signal in Greenland, since 
many boundary conditions responsible for dust defla- 
tion during different climate periods like soil surface 
conditions are difficult to constrain.) If changes in 
dust transport to Greenland can be attributed to ei- 
ther changes in temperature gradient or changes in 
global average temperatures, we might be able to in- 
fer changes in the latitudinal temperature gradient 
from the observed Greenland dust records. This may 
give an additional indication whether the temperature 
gradient information we can obtain from the limited 
amount of ice records in the Northern Hemisphere 
(NH) is sufficient to reconstruct the hemispherical 
temperature conditions. This investigation is specif- 
ically of interest for time periods where land surface 
conditions, which could potentially affect dust defla- 
tion, did not change considerably (vegetation changes 
or addition of glacial outwash), since those changes 
were not included in the simulation. 

2. Dust As Tracer in the GISS GCM 

Dust has been included as a dynamic tracer in the 
GISS atmospheric GCM (4° x 5°horizontal resolu- 
tion, nine vertical layers). This parameterization and 
some results thereof have been described by Tegen 
and Miller [1998]. Dust sources, transport, and de- 
position were computed with a 1-hour time step. In 
the GCM, dust emissions are computed as a func- 
tion of vegetation cover, surface wind speed, and soil 
moisture [Tegen and Fung, 1994], Dust deflation is al- 
lowed in areas labeled by Matthews [1983] as deserts 
or sparsely vegetated regions. The GCM transports 
four size classes of dust, with size ranges of 0.1-1, 1-2, 
2-4, and 4-8 pm as independent tracers. Sizes below 
1 pm were transported as one size class because they 
are not strongly fractionated by gravitational settling. 
Particle sizes larger than 8 pm were not included in 
this calculation, since large particles fall out quickly 
and are not important for long-range dust transport . 


Surface distributions of clay (particles smaller than 
1 pm) and small silt (particle radius between 1 and 
10 pm) were derived from a global soil texture data 
set [Zobler, 1986; Webb et al., 1991], Measurements 
show that above a critical threshold velocity (below 
which no dust deflation takes place), dust fluxes into 
the atmosphere depend on the third power of surface 
wind speed [ Gillette , 1978] which is the most critical 
parameter for calculating dust emissions on the global 
scale. In the model the dust flux in those areas, where 
the surface conditions allow dust deflation, follows 

q a - C(u - u tr )u 2 , ( 1 ) 

where q a is the dust flux from the surface in pg 
m s _l , u is the surface wind speed in m s -1 , and 
u t r * S a threshold velocity. For wind speeds below 
this threshold, no dust deflation takes place. We 
used the dimensional constant of C = 2 pg s 2 nr 5 
for clay particles (< 1 pm) and C = 5 pg s 2 m“ 5 
for silt particles (1-8 pm) to describe dust deflation 
[ Tegen and Miller, 1998]. Because the dust emis- 
sions increase nonlinearly with surface wind speed, 
peak wind speed events are responsible for a major 
part of dust deflation. High-wind events in the GCM 
are less frequent compared to the previously used Eu- 
ropean Centre for Medium-Range Weather Forecasts 
(ECMWF) surface-wind products with a spatial res- 
olution of 1.125° x 1.125° [Tegen and Fung, 1994], 
This leads to an underestimate of dust emissions in 
the GCM experiments. To reduce this difference, we 
chose a threshold velocity for each land grid box in the 
GCM that results in the same dust fluxes compared 
to the off-line model with ECMWF surface winds (for 
which the threshold velocity was 6.5 m s~‘ at all loca- 
tions) [ Tegen and Fung, 1995], The resulting thresh- 
old velocities for GCM surface wind speed vary be- 
tween 4 and 10 m s -1 . We chose to vary the thresh- 
old velocities at each gridbox instead of varying the 
emission factor, because this way the number of dust 
events occurring per year at each gridbox would be 
similar to the off-line results that used the ECMWF 
surface wind product. 

Dust is removed from the atmosphere bv gravita- 
tional settling using Stokes law (size-dependent set- 
tling velocities), turbulent mixing in the first model 
layer, and subcloud washout calculated using GCM 
precipitation. A detailed description of this param- 
eterization and the validation of the results under 
present-day conditions can be found in the work of 
Tegen and Miller [1998]. For these experiments, only 
Asian dust sources are included; that is, dust fluxes 
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pared with the standard deviation than the increase in 
dust concentrations by the increase in the SST gradi- 
ent. Figure 3 shows for case B similar concentrations 
near the source areas compared to the control case 
0, but the concentrations over North America of dust 
transported eastward from Asia are smaller; that is, 
in this case, less dust is transported away from the 
source region. The response of dust concentrations 
to changes in the SST gradient is therefore not sym- 
metrical, the sign of the change in the gradient does 
influence the strength in the dust response. 

An additional decrease in global average SST to- 
gether with the increased temperature gradient (ex- 
periment C) further increases the dust concentration 
at Greenland. However, the increase in the SST gradi- 
ent (experiment A) increases the overall dust concen- 
trations stronger compared to the control experiment 
than a decrease in global average SSTs alone. This 
becomes clear when considering the effect of SST de- 
crease alone by subtracting the results of experiment 
C from experiment A (Figure 3, top middle panel). 
For this case, the concentration of dust above Asia 
is of the same order of magnitude as for the control 
case (experiment 0), but the dust is transported far- 
ther across North America. For experiment C the 
combination of these effects results in dust concen- 
trations near the Asian source areas being similar to 
case A but also in slightly increased dust transport 
and concentration over North America and Green- 
land. There, the dust concentration in experiment 
C is « 20% higher than for experiment A (see Table 
1). This additional increase in Greenland dust is less 
significant than the increase of dust in experiment A 
compared to the control case, considering the large 
standard deviation. 

For the opposite case, if the global mean temper- 
atures are increased additionally to a decrease in lat- 
itudinal temperature gradient (experiment D), the 
dust concentration at Greenland is further reduced 
by about 30% compared to experiment B (Table 1). 
If only the effect of changes in global SST is con- 
sidered compared to the change in the SST gradient 
only (by subtracting the results from experiment B 
from experiment D, middle panel in Figure 3), the 
dust concentrations over both Asian source areas and 
dust transport across North America are slightly re- 
duced compared to the control experiment. For case 
D this change, together with the change due to the 
decrease in SST gradient, leads to a noticeable reduc- 
tion of dust concentrations over the source areas and 
the dust transported eastward. 


Because soil-surface conditions (e.g., changes in 
vegetation cover or land ice) did not change in these 
calculations, dust concentration changes can be caused 
by either changes in dust source strengths or changes 
in transport. Such changes, in turn, can be caused by 
changes in surface wind speed, changes in wind direc- 
tion, or changes in precipitation, which can impact 
both washout rate of the dust during its transport, 
and soil moisture in the source region, changing the 
source strength of dust. 

3.3. Changes in Dust Source Strengths 

Table 1 also summarizes the dust source strengths 
for the different GCM experiments. For increased 
SST gradients the source strength for eastern Asian 
(China) sources is increased by a factor of 3 com- 
pared to the control experiment, while the central 
Asian source is increased by a factor of % 2. This 
corresponds well to the factor 2-3 increase in dust 
concentrations over Greenland for this case. On the 
other hand, a decrease in the SST gradient (case B) 
does not lead to a significant decrease in dust source 
strengths in Asia, as would have been expected if the 
response of the dust to the SST gradient change had 
been symmetrical. 

For experiment C the eastern and central Asian 
source strengths of dust do not increase with an addi- 
tional decrease in global mean SSTs, in fact, they are 
slightly lower than for experiment A. This indicates 
that for an increased gradient the increase in dusti- 
ness at Greenland is controlled by changes in the dust 
source strength, assuming the transport path is the 
same as in the control experiment. With additional 
colder temperatures the additional increase in Green- 
land dust must be caused by changes in transport or 
dust deposition. For the opposite case, an increase in 
global SSTs in addition to a decreased SST gradient 
(case D) leads to an ^ 30% decrease in central Asian 
dust sources compared to case B (increased SST gra- 
dient only), while the Chinese dust sources remain ef- 
fectively unchanged. This decrease is small compared 
to the standard deviation, however. 

This change in dust source strength can be caused 
by changes in surface wind speed, which influences 
the dust deflation, or by changes in precipitation and 
soil moisture, which determines whether a given grid- 
box can act as a dust source (if the vegetation cover 
allows for dust deflation). Plates la-i show differ- 
ence maps for the dust source strengths, surface wind 
speed, and precipitation for experiment A minus ex- 
periment B (the difference between the experiments 
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dust concentrations found in Greenland for the GCM 
integrations with changed average SSTs compared to 
changed latitudinal temperature gradients (Table 1), 
while differences in the dust source strengths cannot 
explain this difference (comparison of experiments A 
and C). 

Dust transport may be affected not only by changes 
in the removal rates but also by changes in trans- 
port pathways. Northward dust transports for exper- 
iments A - D and control experiment 0 are shown in 
Figure 4 for the zonal mean. To compare the differ- 
ences in the transport without including the changes 
in dust source strength and removal rates, the val- 
ues for northward dust transport were divided by the 
total dust content of the atmosphere for each exper- 
iment. The northward transport of dust does not 
considerably change for the different SST boundary 
conditions. For experiment A (increased SST gra- 
dient) the northward dust transport between % 40° 
and 70°N increases compared to the control exper- 
iment, for experiment C (additional cooling) there 
is only a small increase in northward dust transport 
compared to the control experiment. The changes 
in Greenland dust concentrations in the other experi- 
ments can therefore be explained by changes in source 
strengths due to changes in zonal wind speed and by 
changes in washout rates due to changes in precipita- 
tion strengths alone. For the increased temperature 
gradient (without changes in the global mean SST) an 
increase in northward transport also may be a factor 
causing the increased dust in Greenland. However, 
since for this case the factor of increase in Green- 
land dust concentrations agrees well with the factor 
of increase in dust source strength (see Table 1), we 
conclude that the changes in dust transport pathways 
are of comparatively minor importance. 

3.5. Changes in Dust Deposition 

The deposition of dust at Greenland is of interest, 
since the dust signal in ice cores reflects changes in 
the deposited dust on the ground rather than airborne 
dust concentrations at the ice core site. However, as 
mentioned above, the model underestimates the pre- 
cipitation at Greenland and therefore overestimates 
dust concentrations in precipitation and may underes- 
timate dust deposition fluxes at this location. Table 2 
shows total and wet dust deposition at the GISP site 
for the GCM experiments together with the annual 
precipitation and annual mean dust concentration in 
precipitation (based on monthly averages). The wet 
deposition at this site is only about 50% of the total 


(wet plus dry) deposition for cases A and C, 60% for 
the control experiment, and « 70% for cases B and 
D. This reflects the differences in precipitation in the 
model, which is by a factor of 2 smaller for the colder 
climate (experiment C), and by a factor of 2 higher 
for the warmer climate (experiment D) compared to 
the control experiment: for higher precipitation rates 
the wet deposition dominates the dry deposition at 
this remote location. Generally, dust deposition for 
the different GCM experiments follows the trend for 
dust concentrations, where the increase in SST gra- 
dient (experiment A) leads to an increase by a factor 
of 2-2.5 in the dust deposition (and concentration 
in precipitation), while the deposition flux for exper- 
iment C is not significantly higher than for case A. 
On the other hand, the dust concentration in pre- 
cipitation is 30% higher for case C, which is due to 
the differences in the precipitation. The experiments 
with the decreased SST gradient also show a similar 
trend for deposition fluxes as the dust concentrations 
at Greenland, with the total deposition fluxes being 
about a factor of 2 smaller compared to the control 
run. These results indicate that even with the differ- 
ences in precipitation at the Greenland location, the 
changes in deposition fluxes are similar to the changes 
in dust concentrations at this location. 

3.6. Spring Conditions 

As mentioned above, the dust signal at Greenland 
in the GCM shows a summer maximum, although the 
Asian dust source has a spring maximum. This dis- 
crepancy may be due to a high-latitude source that is 
active in the NH summer. Even though Asian sources 
north of 50°N contribute less than 10% and sources 
north of 60°N contribute less than 1% to the total 
Asian dust emission, such sources may disproportion- 
ally influence dust concentrations in Greenland. No 
measurements exist to prove or disprove the existence 
of a dust source at those high-latitude locations. Be- 
cause of lack of measurements the existence of such a 
dust source cannot be excluded. 

To evaluate whether the results presented here 
are still valid for the case that the model did incor- 
rectly predict the existence of such a high-latitude 
Asian source, we investigated the GCM results for 
NH spring additionally to the annual averages. Dur- 
ing NH spring the high-latitude Siberian source is not 
active, while dust emissions in Asia are largest. Table 
3 summarizes the results for the Asian dust sources 
and Greenland dust concentration for NH spring. For 
case A, the increase in springtime dust concentration 
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tributed to greater dust concentrations, although gra- 
dient changes of the values given here were not suf- 
ficient to produce the order of magnitude increase 
found in the GISP ice core. Other obvious differences 
that were not taken into account in these calculations 
could explain these differences: changes in surface 
conditions (glacial outw r ash and vegetation decrease) 
would be expected to fundamentally increase the dust 
source strengths in Asia and subsequent transport to 
Greenland. Furthermore, the existence of large land 
ice sheets in the Northern Hemisphere could well have 
altered circulation patterns. Only a specific simula- 
tion with these boundary condition changes can act 
to provide a full assessment of the magnitude of the 
latitudinal gradient change applicable to the LGM. 
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Figure 1. Seasonal mixing ratios (in fig dust/kg air) of Asian dust aerosol for the control exper- 
iment 0 for the first model layer. 

Figure 2. Seasonality of dust at Greenland for the control experiment 0 for deposition (in mg 
m -2 d* 1 ), first layer concentration (in //g/m 3 ), dust concentration in precipitation (in mg/kg 
water), and precipitation (in mm d _l ). 

Figure 3. Annually and vertically averaged dust mixing ratios for experiments A-D and con- 
trol experiment 0, together with dust concentrations for only decreased (C-A+0) or increased 
(D-B+E) global SSTs (in fig dust/kg air). 

Figure 4. Zonal mean northward dust transports for experiments A-D and control experiment 
0, relative to the total dust content of the atmosphere. 

Plate 1. Differences in dust sources (a,b,c), zonal wind speeds (d,e,f), and precipitation (g,h,i) 
for changes in the latitudinal SST gradient (experiment A minus experiment B) (a,d,g), colder 
conditions without changes in the temperature gradient ( experiment C minus experiment A) 
(b,e,h), and warmer conditions without changes in the temperature gradient (experiment D mi- 
nus experiment B) (c,f,i). 


Figure 1 . Seasonal mixing ratios (in fig dust/kg air) 
of Asian dust aerosol for the control experiment 0 for 
the first model layer. 

Figure 2. Seasonality of dust at Greenland for the 
control experiment 0 for deposition (in mg m -2 d _1 ), 
first layer concentration (in fi g/m 3 ), dust concentration 
in precipitation (in mg/kg water), and precipitation (in 
mm d _1 ). 

Figure 3. Annually and vertically averaged dust mix- 
ing ratios for experiments A-D and control experiment 
0, together with dust concentrations for only decreased 
(C-A+0) or increased (D-B+E) global SSTs (in fig 
dust/kg air). 

Figure 4. Zonal mean northward dust transports for 
experiments A-D and control experiment 0, relative to 
the total dust content of the atmosphere. 

Plate 1 . Differences in dust sources (a,b,c), zonal wind 
speeds (d,e,f), and precipitation (g,h,i) for changes in 
the latitudinal SST gradient (experiment A minus ex- 
periment B) (a,,d,g), colder conditions without changes 
in the temperature gradient ( experiment C minus ex- 
periment A) (b,e,h), and warmer conditions without 
changes in the temperature gradient (experiment D mi- 
nus experiment B) (c,f,i). 


Table 2. GCM Results for Annual Average Precipitation and Dust Deposition at Greenland 



Depos. Flux 
mg/m 2 /yr 

Depos. Flux (Wet) 
mg/m 2 /yr 

Precipitation 

rnm/yr 

Cone. (Precip.) 
Hg/kg Water 

C (A ± colder) 

52 ± 10 

22 ±5 

36 ±9 

490 ± no 

A (increased gradient) 

51 ± 12 

23 ± 7 

57 ± 16 

380 ± 90 

0 (control) 

22 ±6 

14 ± 4 

83 ±22 

160 ± 58 

B (decreased gradient) 

12 ± 3 

8 ± 3 

98 ±23 

79 ± 22 

D (B ± warmer) 

10 ±3 

7 ± 3 

150 ± 30 

53 ± 18 
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Table 3. GCM Results for Dust Sources and Dust Concentrations at Greenland for NH Spring 



China Source 
mg/m 2 /yr 

Central Asia Source 
mg /m 2 /Seas 

Cone. (Layer 1) 
M g/'n 3 

Avg Mixing Ratio 
Mg/kg Air 

C (A + colder) 

5.3 ±2.6 

20 ±8 

0.71 ±0.11 

1.6 ±0.2 

A (increased gradient) 

13. ± 5. 

31 ±10 

1.1 ±0.24 

2.2 ± 0.3 

0 (control) 

3.7 ±2. 

15 ±6 

0.27 ± 0.09 

0.9 ± 0.2 

B (decreased gradient) 

3. ± 1.8 

8 ± 4 

0.16 ± 0.04 

0.6 ± 0.1 

D (B 4- warmer) 

1.8 ±1 

7 ± 3 

0.10 ±0.03 

0.5 ±0.1 
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